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Abstract 

In view of current interest in relativistic spin-one systems and the recent work on the 
Dirac Oscillator, we introduce the Duffin-Kemmer-Petiau (DKP ) equation obtained by using 
an external potential linear in r. Since, in the non-relativistic limit, the spin 1 representation 
leads to a harmonic oscillator with a spin-orbit coupling of the Thomas form, we call the 
equation the DKP oscillator. This oscillator is a relativistic generalisation of the quantum 
harmonic oscillator for scalar and vector bosons. We show that it conserves total angular 
momentum and that it is exactly solvable. We calculate and discuss the eigenspectrum of 
the DKP oscillator in the spin 1 representation. 


1 The DKP Oscillator 

The focus of attention in this paper is to generalise the concept of the quantum harmonic 
oscillator to relativistic vector bosons. 

For a free scalar or vector boson of mass m, the relativistic DKP equation [1] is 

(c/3 ■ p + mc 2 )if = ih/3 0 -^ (1-1) 

where the internal variables ^ (/r = 0,1, 2, 3) satisfy the commutation relation 

pL , f3 l '(3 x + p x p v p lt = g tlu P x + g uX P tl - ( 1 * 2 ) 

In the spin 1 representation, the are 10 x 10 matrices while the dynamical state 4> is a 
ten-component spinor. 

For the external potential which we introduce with the non-minimal substitution 

p — > p — imujr] 0 r, (1-3) 

where uj is the oscillator frequency and r/° = 2/?° 2 — 1, the DKP equation for the system is 

[c/3 ■ (p — imijjifr) + mc 2 ]ift = (1-4) 

This external potential, which is of Lorentz tensor type, does not conserve the orbital and spin 
angular momenta, since 

[/3t/° • r, L] = -i(f3rj° A r) and [prf • r, S] = i(/3r/° A r), (1.5) 
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but it does conserve the total angular momentum J = L + S. 

In the spin 1 representation of eq.(4), the dynamical state ip is chosen as the 10-component 
spinor 

(CA 

C 2 (1.6) 
\Cs) 

so that, for stationary states, the equation of motion eq.(4) decomposes into 

{ me 2 (f = icp~ * B 

me 2 A — EB — cp + AC , , 

mc 2 B = EA + icp+<p 
mc 2 C — — cp” A A 

where p ± = p ± imwr. Since A is the 3-component spinor analogous to the Dirac upper 
component, we seek the wave equation for A. It is straightforward to eliminate < p , B and C in 
favor of A so that one gets 

(E 2 — m 2 c 4 )A — [c 2 (p 2 + m 2 w 2 r 2 )- 3^mc 2 - 2^mc 2 L*s]A ^P + {P _, [p + A(p“AA)]} (1.8) 

where L is the orbital angular momentum and s the 3x3 spin one operator. In the non- 
relativistic limit e <C me 2 , the fourth term in eq.(1.8) becomes negligible, since it is of order 
1/m 3 , so that the wave equation for A can be written 


V>(r) - 


< "f \ 

A(r) 

B(r) 

VC(r)/ 


with A 



R = 



p2 i 3 

e A ~ f- (- -mw 2 r 2 tux! — %ujh ■ si A 

2m 2 2 J 


(1.9) 


which characterises the usual harmonic oscillator in addition to a spin-orbit coupling, absent 
for scalar DKP bosons, of strength —hu?. Note that the strength of this coupling is half the one 
obtained from the Dirac oscillator [2]. 

Since the spin 1 representation of eq.(1.4) leads to the usual three-dimensional (3D) oscillator, 
in the non-relativistic limit, we refer to the system it describes as the Duffin-Kemmer-Petiau 
oscillator. 


2 Solution to the vector DKP oscillator problem 


For the S = 1 central field problem, the general eigenfunction we use takes the form [3] 


^jm{ r) 


r Z L G nJL (r)Y ft, (SI) 


( 2 . 1 ) 


Putting V’JM i 11 * 0 eq.(1.4) results in ten coupled radial differential equations which can be 
decoupled into two sets associated with ( — 1) J and (— 1)' /+1 parities. We call the (— 1) J solutions 
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natural-parity (or magnetic-like) states while we refer to the ( — 1) J+1 solutions as unnatural- 
parity (or electric-like) states. With the notation 


the set associated with ( — l ) 7 parity is 


, d J + 1 mujr 

%c[ Tr ~ + ~r 


- (j 


he 


d J + 1 inter 
dr r h 


d J inter 
dr + r ^ % 


d_ 

dr 


) = R± i, 

R = F,G,H 

( 2 . 2 ) 

m.c 2 G 0 


(2.3) 

) F ° = 

--—mc 2 H\ 

G 

(2.4) 

) F a = 

~—m.c 2 H_ x 

(2.5) 

J 

<*.) 

J rater 

\ 1 / , X 

( 2 . 6 ) 

r h 

J H_\ = — (me 2 F a - EGo) . 


For unnatural parity states, the radial differential equations are coupled in the following 
way : 


he 


d .7 + 1 nuor 
dr r fi 


. d J niter 

%C fc + r fT 


= ~^{m.c 2 F x -EG x ) 
j H a = -+ (rnc 2 F„ x - EG_ X ) 


Id .7 + 1 j liter \ „ id J niter _ 1 9 

G 7 - + + — T~ Fi - ctj + — — F_, = — m,c 2 Ho 

dr r Ti I dr r Ti he 


(2.7) 

(2.8) 

(2.9) 


Id J + 1 rmer \ 

he — <f> 

dr r hi 


, . d, J inter . . 

he 1 — + — <f> 

dr r h 


— — (mc 2 G x - EF X ) 

nj v 7 

-I (mc 2 G_i - EF_ X ) 

G v 7 


( 2 . 10 ) 

( 2 . 11 ) 


( d J + 1 mojr \ . / d J mu>r \ 1 2 

7- + + -7— Gi+oh + -7— G-1 = r-mc 2 G (2.12) 

dr r h I \drrnl he 

To obtain the exact solution for the magnetic-like states, we eliminate G 0 , H± in favor of F 0 
in eq.(2.6). This yields the eigenvalues [4] 


2 me 2 


(Efj j - m 2 c 4 ) = (N + 1 )hte 


(2.13) 
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with the principal quantum number N = 2n + J (n is the radial quantum number). Note that 
the oscillator levels are equidistant and degenerate; the zero-point energy differs here from the 
one we found for the scalar DKP bosons. 

The exact eigenvalues of the radial equations associated with unnatural-parity states can be 
shown [4] to be 

2 h iE * ~ mV > = {N + i>*“ + 7(7 + A (2 - 14) 


where 



1 + 


Oi tux) 

a 0 me 2 




1 

2 


(2.15) 


with a 0 = (2 J + l) 2 , a\ = 4 J(J + l)(2iV + 3) and a 2 = 4 J 2 (J + l) 2 where N, a positive integer, 
is the principal quantum number. 

As shown in eq.(23), the energy of the DKP oscillator in unnatural parity states involves the 
usual 3-dimensional harmonic oscillator energy, a second term proportional to J(J + 1) which 
appears as some kind of rotational energy and a third energy contribution A which is a compli- 
cated function of the oscillator frequency, J and N with no obvious physical interpretation. 

In the limit where the oscillator frequencies are such that hu> <C me 2 , keeping only the 
first-order term in u> in eqs.(2.14-15) leads to 


A^(Ej-mV) = e+. ~(N-J+l)hu, (2.16) 

= ~(N + J + 2)hu (2.17) 

2 me 2 

This is best illustrated in fig.l which shows, for fixed values of N and J, the variations of the 
relativistic and non-relativistic eigen-energies with hcj/mc 2 . 



hco/mc 

Figure 1: Variation of the DKP and non-relativistic oscillator energies with LJ. 


This shows that our solutions have the correct non-relativistic limits since the levels in 
eqs.(2.16-17) are those of a usual 3D non-relativistic oscillator with a spin-orbit coupling of 
strength -hu>. In this limit, they could have also been obtained directly from eq.(1.9). Fur- 
thermore, taking this limit suggests the interpretation of the E + and E _ energies as “spin-orbit 
partners”, E+ being associated with J = L + 1 and E_ with J = L — 1. 
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The unnatural parity E + levels for N < 9 are presented in fig.2 alongside the e+ r and the 
( N + |)hu; levels for reference. For a given value of N - J, all the non-relativistic energy levels 
( N,J *■) (with N — J odd and J < N) are infinitely degenerate. This accidental degeneracy is 
not present in the exact DKP oscillator E + eigenspectrum whose levels are found to cluster in 
several groups of states belonging to the same N — J oscillator shell. 



Figure 2; DKP and non-relativistic spectra associated with J = Ij t 1 for N < 9. The dotted lines between the DKP 
and non-relativistic oscillator levels link states with the same quantum numbers ( N , «/")• 


The E- eigenspectrum is now presented in fig. 3 together with the non-relativistic e n r energy 
levels for N < 9. 



While all the non-relativisitic (N,J*) levels associated with the same N + J oscillator shell 
are degenerate, with a finite degeneracy in this case, their relativistic analogues are not. The 
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exact DKP oscillator states are found to cluster into bands of states belonging to specific values 
of N + J. 

For a more quantitative analysis of the E- bands of the DKP oscillator spectrum, we plotted 
in fig.4 the E _ energy levels, belonging to the N + J = 49 band for instance, against J(J + 1) 
for different oscillator frequencies. 



J(J+1) 

Figure 4: Energy levels of the N + J = 49 band as a function of J(J p 1 ) for different oscillator frequencies. 

It is indeed remarkable that the DKP oscillator energies constitute nearly perfect rotational 
bands. There are deviations from the rotational patterns at low angular momenta. These single 
particle rotational bands are of the finite type since for N P J fixed they terminate at some 
Jmax ■ The effective rotational moments of inertia are sensitive to the oscillator frequencies since 
the slopes of the bands are found to vary substantially with increasing uj. 



J ( J + 1 ) 

Figure 5: Energy levels of the N p J — 23, 29, 35, 41, 49 band as a function of t/(<7 + 1) for h(J = 0.2GeV\ 
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Fig.5 alternately represents the energies of five different N + J bands for a specific oscillator 
frequency as a function of J(J + 1). The DKP oscillator energies now lie on rotational bands 
whose slopes hardly change with N + J which implies that the effective rotational moments of 
inertia are rigid and insensitive to such variations. 

Of course, it should be pointed out that these rotational bands are unlike the usual ones 
where the levels are associated with the same intrinsic motion but different angular momenta. 
Here the single particle states involve different radial as well as rotational motions. Note that 
this behaviour is not particularly tied to this DKP oscillator. Buck [5] also found that, when 
solving the Schrodinger equation for deep, bell-shaped potentials, the levels with a fixed value 
of N = 2ti + £ ( these states are degenerate for a harmonic oscillator ) lie on a straight fine when 
plotted against £(£ + 1). Geometric arguments in terms of the shapes of the potentials which can 
give rise to these rotational-like bands have been put forward to explain this behaviour [5] [6]. 


3 Conclusion 

We have introduced a new potential in the DKP equation. Since, in the non-relativistic limit, 
the DKP equation of motion leads to the usual harmonic oscillator with a spin-orbit coupling 
of the Thomas form, we call the system a DKP oscillator. This oscillator is a relativistic 
generalisation of the quantum harmonic oscillator for vector bosons. We have shown that it 
conserves the total angular momentum, that it is exactly soluble and we have computed and 
discussed its eigen-solutions. 

The renewed interest in the Dirac oscillator has generated studies of its group theoretical 
properties [7] and hidden supersymmetry [8] [9] among others. Such investigations of the DKP 
oscillator would be most useful to gain further insight into the physical meaning of this oscillator. 

This study is on the other hand relevant to the work on relativistic equations for two fermions 
and particularly to those of Krolikowski’s type [10]. Since they tend to the DKP equation in 
the point-like limit of tightly bound-states, exact solutions of the latter may provide usefid 
information about this class of relativistic two-body equations. 
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